Vortex sorter for Bose-Einstein condensates 
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We have designed interferometers that sort Bose-Einstein condensates into their vortex compo- 
nents. The Bose-Einstein condensates in the two arms of the interferometer are rotated with respect 
to each other through fixed angles; different vortex components then exit the interferometer in dif- 
ferent directions. The method we use to rotate the Bose-Einstein condensates involves asymmetric 
phase imprinting and is itself new. We have modelled rotation through fixed angles and sorting 
into vortex components with even and odd values of the topological charge of 2-dimensional Bose- 
Einstein condensates in a number of states (pure or superposition vortex states for different values 
of the scattering length). Our scheme may have applications for quantum information processing. 

PACS numbers: 03.75.-b,03.75.Kk,05.30.Jp 
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Introduction. One of the central characteristics of a su- 
pcrfluid such as a Bose-Einstein condensate (BEC) is the 
presence of quantized vortices. Vortices have been gen- 
erated experimentally 0, 0] by stirring the BEC, very 
similar to the rotating-bucket experiment in Helium [3|. 
The detection of vortices in BECs is typically made by 
a direct observation of the vortex core or by interference 
experiments 4J. In present experiments the charge, m, of 
the vortex is in principle known as the initial rotation fre- 
quency which stirs the cloud is known. In this paper we 
show how to sort vortices when the charge is not known. 

Vortices have attracted considerable interest both ex- 
perimentally and theoretically, mainly because of their 
inherent many-body character and the connection to fluid 
dynamics. In addition, optical vortices in single pho- 
tons have recently been used to carry information - and 
in particular quantum information |5|. Light is an ex- 
cellent carrier of information over large distances as the 
photons travel very fast and do not easily interact with 
each other. For the same reason, photons are not very 
well suited for storing the information for longer times. 
This is where atoms would be better suited as a medium 
for storing information, especially quantum information. 
Optical vortices, special cases of light with orbital angu- 
lar momentum, can carry huge amounts of information as 
there is in principle no limit to the quantized angular mo- 
mentum acting as the information carrier. Transferring 
this angular momentum to atoms would constitute a way 
to store the information |(|. If such a storage device is 
to work with atoms we need a way to manipulate atomic 
states, and in particular vortex states, in an efficient and 
useful way. It is therefore important to know the mecha- 
nisms behind the vortex dynamics and more importantly 
how to manipulate the vortex states in order to be able 
to make any kind of readouts from the trapped quantum 
gas. In this paper we study theoretically the application 
to BECs of ideas borrowed from optical vortex sorting • 
Vortex sorter. If a vortex (in light or in a BEC) of charge 
m = 1 is rotated through 180° about its centre, it changes 



phase by 7r (and, in the simplest case, is unchanged in 
any other respect). If, on the other hand, a vortex of 
charge m = 2 is rotated through 180°, its phase is un- 
changed. The two cases discussed above are in fact rep- 
resentative for all vortices with odd and even charges, 
respectively. This effect has been used in an optical two- 
arm interferometer which rotates the beams in the two 
arms with respect to each other to route vortices accord- 
ing to their charge into one of the interferometer's two 
exit ports Q. When the beams are re-combined, even- 
charge vortex components interfere constructively in one 
interferometer port and therefore exit the interferometer 
through it, while odd-charge vortex components interfere 
destructively in that port and therefore exit the interfer- 
ometer through another port (in which even-charge vor- 
tices interfere destructively). The vortices exiting from 
the two ports can be sorted further in similar interferom- 
eters, but with different relative rotation angles For 
example, vortices with even charges (m = 0,±2,±4, ...) 
can be sorted into those whose charges are respectively 
integer and half-integer multiples of 4 (m mod 4 = 
or m mod 4 = 2, respectively). In some cases, uniform 
phase offsets in one arm are required ||. 

By using Bragg pulses, it is possible to coherently split 
and re-combine a BEC as would be required within a 
two- arm interferometer (figure ^) 9], and using specially 
designed light pulses a BEC could be rotated through any 
given angle (see below). We numerically examine here a 
vortex sorter created by combing these two elements as 
shown in figure ^ 

Specifically, we model a 2-dimensional BEC that is 
split into two identical BECs, which are then rotated 
with respect to each other through 180°, and finally su- 
perposed. The wave function of the original BEC, 'F, is 
split according to 
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The two BECs are rotated with respect to each other 
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FIG. 1: Schematic of a vortex-sorting Bragg-pulse interfer- 
ometer. A stationary input BEC is first split by a n/2 Bragg 
pulse (left; the horizontal axis is time). The two resulting 
BECs are in different momentum states and move apart. A 
7r Bragg pulse swaps the momentum states, so that the two 
BECs move together again. We refer to the two different tra- 
jectories in this space-time diagram as two 'arms' of the inter- 
ferometer. The BECs in the two arms are flipped vertically 
and horizontally, respectively, which corresponds to a rela- 
tive rotation through 180°. A second ir/2 Bragg pulse mixes 
the two BECs such that even and odd vortex components 
have different momentum states - they exit the interferome- 
ter through different 'ports'. 



through mirroring the two BECs with respect to the x 
and y axis, respectively. As explained below, this can 
be achieved through imprinting specific phases onto the 
BECs at specific times. In combination with the time 
evolution between the phase-imprinting events, this re- 
sults in wave functions and ^' 2 . We simulate this 
time evolution according to the time-dependent Gross- 
Pitaevskii equation which we write in the form 



— V 2 + -muj 2 r 2 + g\ip\ 2 
2m 2 
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where m is the mass of each atom, ui is the trap frequency 
and g is the non-linear coefficient, which is given by 
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Here d is the effective thickness of the BEC in z direc- 
tion and N is the number of atoms. The motion can be 
considered two-dimensional if the chemical potential fx 
of the trapped cloud is smaller than the corresponding 
trapping energy huj z in the z direction. Finally, the two 
wave functions "F^ and are superposed according to 
the equations 
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This model can represent various interferometers, all 
of which are idealised in some respects. For example, a 
Bragg-pulse interferometer with rotation in the arms (fig- 
ure^), is idealised as follows. Firstly, the Bragg pulses are 



assumed to be perfect, that is acting according to equa- 
tions and which describe perfect ir/2 pulses with 
the exception that the two states <Fi and have differ- 
ent momenta. The n pulses, which are also required in 
the Bragg-pulse-interferometer scheme, and which swap 
the BECs between the two states, are also assumed to be 
perfect (which is consistent with experiments in which 
fringe visibilities close to 1 were achieved in Bragg-pulse 
interferometers [13 )• Secondly, the interaction between 
the BECs in the different arms is neglected. To the best 
of our knowledge, no experimentally realisable situation 
is exactly represented by this, but some are represented 
better than others. A Bragg-pulse direction that sepa- 
rates the planes of the two BECs, for example, should 
lead to less interaction between the BECs than a Bragg- 
pulse direction that move the two planes across each 
other; however, applying the light pulses for rotation to 
the two arms separately is potentially difficult in this ge- 
ometry. Thirdly, the arm length is just that required for 
rotation; we have made no allowance for any additional 
time it might take for the BECs to separate sufficiently 
such that they can be rotated independently and sub- 
sequently recombined. However, in analogy to optical 
imaging, lens light pulses might be able to return the 
BEC into an earlier state, thereby effectively shortening 
the arms. Other interferometer types that approximate 
our idealised model include, for example, those that split 
a BEC into two by putting them into different internal 
states and manipulate the two BECs independently 
through phase imprinting with light pulses with different 
detunings. 

Figure [5] shows examples of sorting a BEC into its 
'even' and 'odd' vortex components, which is demon- 
strated for pure vortex states as well as superpositions. 
These simulations - indeed all the simulations in this pa- 
per - were performed over an area of 14 x 14 (in units 
of yjTi/muj) on a 256 x 256 grid of wave function ampli- 
tudes. Figure |3] shows the fraction of the original BEC 
that exits the interferometer in the correct port - a mea- 
sure of the quality of the sorting - as a function of the 
non-linear coefficient g. It can be seen that the scheme 
works better for small values of g. 

Rotation of BECs. Several methods already exist for ro- 
tating a BEC through a given angle, using, for example, 
an external magnetic field [l4|, a careful arrangement of 
laser beams 0, or Bragg pulses 0, 0. We describe 
here a novel method based on an optical analogy [l7| . 

Our method is based on the fact that mirroring at one 
axis (or plane in 3 dimensions), followed by mirroring at 
another axis, which is rotated with respect to the first 
axis by an angle a, is equivalent to a rotation through 
an angle 2a about the intersection between the two axes. 
In analogy to mirroring of a light beam, which can be 
achieved with a pair of identical cylindrical lenses par- 
allel to the mirror axis, each of focal length /, which 
are separated by 2/ (such a configuration is called a it 
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FIG. 2: Simulated vortex sorting in an idealised interferom- 
eter. The columns show the probability density and phase 
(inset at reduced size) of the input BEC (left) and the BEC 
in the 'even' (centre) and 'odd' (right) output ports. The top 
three cases show sorting of pure vortices in m = 1, m = 2, 
and m = 3 states; depending on the value of m as it enters the 
interferometer, the BEC comes either out of the even or odd 
output port. This sorting is not perfect: a small fraction can 
be seen to come out of the 'wrong' output port. The fourth 
case shows sorting of a superposition of vortices with charges 
m = 1 and m = 2, which get split into its vortex compo- 
nents. This figure was calculated for g — 0. Note that most 
of the structure near the edge of most phase plots in this pa- 
per, where the probability density is very low, is a numerical 
artefact. 



mode converter |ljfl), a BEC can be mirrored by a pair 
of correctly spaced cylindrical-lens pulses. These are far 
off-resonant light pulses with a transverse intensity dis- 
tribution that is proportional to the thickness of the cor- 
responding optical cylindrical lenses, that is the intensity 
falls off quadratically in one direction and is constant in 
the other. The effect of each cylindrical-lens pulse is a 
phase change proportional to the local intensity 0, 0] : 
the cylindrical-lenspulses act like phase holograms of 
cylindrical lenses [13j. The phase change due to each 
lens pulse is r 2 /(Atf), where r is the distance from the 
axis of the cylindrical- lens pulse and tf is its focal time 
(the equivalent of the focal length in optical lenses), fn 
this paper we use tf = 0.03 (in units of l/u>), which is 
one of the smallest focal times that satisfies the Nyquist 
criterion for our model, and a time of f<j = 0.06 between 
the lens pulses. Figure 0] illustrates modelled examples 
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FIG. 3: Fraction of the BEC in the correct port as a function 
of the non-linear coefficient, g, and for different values of the 
topological charge. Inset are the density cross-sections in the 
even (left) and odd (right) ports corresponding to some of the 
data points (m = 1, g = 500 and m = 4, g = 250 and 750). 

of rotation of BECs through 180°. 

This scheme does not work perfectly, not even in op- 
tics: a light beam (and, by analogy, a BEC with g = 0) is 
mirrored perfectly only in the limit of cylindrical lenses 
with infinitely short focal lengths 01 ■ Obviously, this is 
not possible, and the result is imperfect mirroring that 
leads to asymmetry and vortex splitting. Another prob- 
lem when using cylindrical-lens pulses to mirror BECs 
with g 7^ is that the BEC can intermittently become fo- 
cussed into a line, which greatly amplifies the non-linear 
effects, which in turn usually lowers the quality of the 
mirroring. 

In the context of the vortex-sorting interferometer, we 
are interested in differential rotation between the BECs 
in the two arms of the Bragg-pulse interferometer. A 
better way of achieving such differential rotation is to 
apply the first two of the four rotation pulses shown in 
figure 0] to the BEC in one arm, and the other two to 
the BEC in the other arm; both BECs are mirrored, but 
with respect to different axes. As demonstrated above, 
we find that - in the spirit of spreading imperfections 
symmetrically and hoping for cancellation - this leads to 
good results for g > 50 for small values of \m\. 
Conclusions. In this paper we have investigated the 
sorting of vortices in BECs using an intcrfcrometric 
technique. Our technique requires the BEC to be ro- 
tated, which we achieve with spatially inhomogeneous 
imprinted phases. If the non-linearity is strong, the effi- 
ciency of the rotation and therefore the efficiency of the 
sorting process are decreased, but with existing experi- 
mental techniques such as exciting Feshbach resonances 
it should be possible to "tune" the non-linear coefficient 
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FIG. 4: Rotating a BEC through 180° with a series of 
cylindrical-lens pulses. Shown at the top are the intensity 
cross sections of four far off-resonant light pulses. The inten- 
sity of the first two pulses falls of quadratically in the y, that 
of the other two pulses falls of quadratically the x direction. 
The pictures below show examples of simulated BECs before 
interaction with the light pulses (left), after interaction with 
the first pair of pulses (centre) , and after interaction with the 
second pair (right), a: pure vortex with m = 1, g = 500. 
It can be seen that the handedness of the phase distribution 
(inset) is reversed after the first pair of lens pulses, and that it 
is back to its original handedness after the second pair. That 
this leads to rotation can be seen more clearly in b, which 
starts off with a superposition of m = 1 and m = 2 vortices 
(g = 0): after interaction with all pulses the probability den- 
sity is indeed rotated through 180°. In both cases, the focal 
time of each cylindrical-lens pulse was tj = 0.03 (in units 
of 1/cj); the pulses passed the BECs at times t = 0, 0.06, 
0.06 (the two middle pulses were merged into one), and 0.12, 
respectively. 

g to a value suitable for reliable vortex sorting. 

The techniques presented here are based on methods 
used in conventional optics. When these methods are 
transferred to BECs, complications arise, but also some 
intriguing new possibilities. In particular non-linearity 
- the origin of most complications - is important when- 
ever information is not only to be stored, but also to be 
processed in computations. 
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